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ABSTRACT 

We determine the one-loop and the two-loop back-reaction corrections 
in the spectrum of the Hawking radiation for the CGHS model of 2d dila- 
ton gravity by evaluating the Bogoliubov coefficients for a massless scalar 
field propagating on the corresponding backgrounds. Since the back- 
reaction can induce a small shift in the position of the classical horizon, 
we find that a positive shift leads to a non-Planckian late-time spectrum, 
while a null or a negative shift leads to a Planckian late-time spectrum 
in the leading-order stationary-point approximation. In the one-loop case 
there are no corrections to the classical Hawking temperature, while in 
the two-loop case the temperature is three times greater than the classi- 
cal value. We argue that these results are consistent with the behaviour 
of the Hawking flux obtained from the operator quantization only for the 
times which are not too late, in accordance with the limits of validity of 
the semiclassical approximation. 
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1. Introduction 



Two-dimensional (2d) dilaton gravity models have turned out to be excellent toy 
models for understanding black hole formation and back-reaction of the Hawking ra- 
diation (for a review see U). The CGHS model has been studied most extensively, 
due to its simplicity. In || it was shown that a thermal Hawking radiation can exist in 
the semiclassical limit, which was done by evaluating the Hawking flux from the trace 
anomaly. This result was later confirmed by a direct calculation of the Bogoliubov 
coefficients which was the 2d analog of the famous Hawking calculation [[|. In a 
further development, explicit solutions which include the back- reaction were found, 
in one-loop || |6], |7] and two-loop approximations. 

In the one-loop case, the trace anomaly method |7|], as well as the point-splitting 
method |§, yield the same Hawking flux as in the zero-loop case, which is a good indi- 
cation that the one-loop back-reaction does not change the Planckian spectrum of the 
radiation. This result is surprising, because one expects that the back-reaction should 
change the zero-loop temperature for a small amount. This is especially puzzling in 
the case of the BPP solution [[7], where the thermal Hawking flux gets turned off 
for late times (although discontinuously), leaving behind a static remnant geometry. 
Another puzzling feature of the BPP solution is that the corresponding semiclassical 
geometry has a horizon which is shifted for a small amount with respect to the clas- 
sical horizon. The horizon shift also appears in the two-loop case, while the operator 
quantization yields an expected result for the Hawking flux, in the sense that the 
zero-loop flux changes under the two- loop back-reaction 0. 

These results make interesting the calculation of the Bogoliubov coefficients, since 
one would like to see what kind of corrections these back-reaction effects would make 
in the spectrum of the Hawking radiation. Another interesting point to be examined is 
that one can evaluate the Hawking flux from the Bogoliubov coefficients and compare 
it to the Hawking flux obtained from the expectation value of the energy-momentum 
tensor operator. Therefore we are going to evaluate the Bogoliubov coefficients for a 
free scalar field propagating on the one-loop geometry of @ and the two-loop geometry 

of i. 

In section two we briefly review the quantization of the CGHS model in the op- 
erator formalism of 0, and define the loop expansion. In section three we evaluate 
the Bogoliubov coefficients for the BPP background geometry and analyse the corre- 
sponding spectrum of the Hawking radiation. The same is done in section five, but in 
the case of the two-loop background geometry of ref . || . In section six we present our 
conclusions. In the appendix we give some theorems which are useful for obtaining 

2 



the asymptotic behaviour of the Bogoliubov coefficients. 



2. Loop-expansion for 2d dilaton gravity 

Our starting point is a classical theory described by the CGHS action 

S = j M <F x ^-g[e^(R + {V<P) 2 + A\ 2 )-\{Vf) 2 ] , (2.1) 

where is a dilaton scalar field, / is a matter scalar field, g, R and V are the 
determinant, the curvature scalar and the covariant derivative respectively, associated 
with a metric g^ on a 2d manifold M. The topology of M is that of R x R. The 
equations of motion can be solved in the conformal gauge ds 2 = —e p dx + dx~ as 



e 



-p 



-<P = _X 2 x + x- -F + -F_ , / /.(./••) • / i.r ) , (2.2) 



where 



F± = a± + &±x ± + J X dy J V dzT±±(z) , (2.3) 
and T±± is the matter energy-momentum tensor 

T±± = kd±fd ± f . (2.4) 

The residual conformal invariance is fixed by a gauge choice p = <f), and the indepen- 
dent integration constants are a + + a_ and b±. 

The quantum theory can be obtained by quantizing the space of classical solutions 
defined by the equations (2.2) and (2.3) 0. This is equivalent to a reduced phase space 
quantization of the action (2.1) 0. In this way one obtains a quantum theory of a free 
massless scalar field / propagating on a flat fictitious background ds 2 = —dx + dx~. 
The dilaton and the conformal factor are given by the expression in (2.2), which is 
considered as an operator in the Heisenberg picture. The matter energy-momentum 
tensor operator is defined as 

T±± = | : d ± fd ± f : , (2.5) 

where the normal ordering in (2.5) is chosen to be with respect to the vacuum for the 
creation and annihilation operators (a k , a^) defined by 



= -^r^r= V^~ lkx± + <^ tkx± ] > (2-6) 

y/2-K Jo \f2uJk 1 J 



where — \k\. 



The physical Hilbert space of the model is the Fock space J-(a k ) constructed from 
a k acting on the vacuum |0). The model is unitary because the dynamics is generated 
by a free-field Hamiltonian 

f°° t 

H— dku k a\a k + E , (2.7) 

J-OO 

which is a Hermitian operator acting on J 7 , where £0 is the vacuum energy. Con- 
sequently the states at t = \{x + + x~) = const, surfaces are related by a unitary 
transformation 

ty(t 2 ) = e-^-^ih) . (2.8) 
The Heisenberg picture is defined as 

tf = e iHt ^{t) , A(t) = e iHt Ae- iHt , (2.9) 

which also serves to relate the covariant quantization expressions to the canonical 
quantization expressions. For example, 

f(t, x) = e iHt f(x)e~ iHt = -L f°° -jj= fa fc e <(te - w * t) + ale'^'^} , (2.10) 

\/2lX J-00 \/2uJ k 

where x = \{x + — x~). 

Given a physical state \l/o; one can associate an effective metric to = e~ tHt ^o 

via 

where e Peff is the effective conformal factor. The geometry which is generated by e peff 
via ds 2 = —e Pef! dx + dx~ makes sense only in the regions of M where the quantum 
fluctuations are small || ||. The condition for this is 

sj\ (e 2 ^) - (e*f I « (e*) , (2.12) 

which defines the limits of validity of the semiclassical approximation. At least per- 
turbatively, the region defined by (2.12) roughly coincides with the weak-coupling 
region @, §. 

The effective conformal factor can be calculated perturbatively by using an expan- 
sion in powers of the energy-momentum tensor, which is equivalent to the expansion 
in matter loops |], || 

00 

((-\ 2 x + x- - F)- 1 ) = e* Y, en * ( 5pn ) > ( 2 - 13 ) 

n=0 
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where F is a c-number function, e *° = — X 2 x + x — F and 5F = F — F . A 
convenient choice for Fo is 

F = <F + ) + (F_> , (2.14) 
since then the lowest order metric is a one-loop semiclassical metric 

e -0o = ^\i x + x - _ _ . (2.15) 

The state \l/o is chosen such that it is as close as possible to the classical matter 
distribution fo(x + ) which creates the black hole. The corresponding classical metric 
is described by 

! n " ^ - X 2 x + A(x + ) - X 2 x + x~ (2.16) 



A 

where 

M(x + ) = A f + dyyTl + {y) , A 2 A = f + dyT° ++ (y) (2.17) 

J — oo J—oo 

and T° + = ^d + f d + f . The geometry is that of a black hole of the mass 



M= lim M(x + ) , (2.18) 

£+— > + oo 



and the horizon is at 



= -A = - lim A{x + ) . (2.19) 
In the limit of a shock-wave matter distribution, for which 

T° + = a5(x + - x+) , (2.20) 

we have 

M(x + ) = \axt9(x + - atf) , A(rr+) = . (2.21) 

A 

The asymptotically flat coordinates (r] + ,ri~) at the past null infinity X~ are given by 

Xx + = e Xr,+ ) x - = -Ae- Xr >~ , (2.22) 
while the asymptotically flat coordinates (a + , a~) at the future null infinity X + satisfy 

\ x + = e Xa+ , X(x~ + A) = -e- A,T_ . (2.23) 

Note that a conformal change of coordinates x ± — > ^ defines a new set of creation 
and annihilation operators (b\,bk) through 
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The old and the new creation and annihilation operators are related by a Bogoliubov 
transformation 

/oo 
dq(b q a qk + blP* qk ) , (2.25) 

and the new vacuum is given by |0^) = S |0). Since initially the geometry should be 
as close as possible to the dilaton vacuum, we take for ty a coherent state 

0-) , (2.26) 

0~^> is the vacuum for the "in" coordinates (2.22), while 

roc 

A= dk\f (k)al k -f*(k)a. k ] , (2.27) 
Jo 

and fo(k) are the Fourier modes of fo(x + ). The choices (2.26) and (2.27) complete 
the initial set up of the quantization method of |J. 



where |0^) = 



3. One-loop spectrum 

A one-loop effective metric can be obtained from the expression (2.15), which was 
calculated in ||. This gives 

e~ p0 = e~^ = C + b±x ± -\ 2 x + x- -^\og\\ 2 x + x-\-\ p dy + (x + -y + ) (Jf^J ■ 

V (3-1) 

The expression (3.1) can be also obtained as a solution of the equations of motion of 
an effective one-loop action 

S eff = S -^J d 2 x^jRU ~ X R-k( d 2 xy/^{R<j) - (V0) 2 ) , (3.2) 

where S is the CGHS action (2.1)0. By choosing C = -\k[1 - log(«/4)] one can 
obtain a consistent semiclassical geometry . In the case of the shock- wave matter 
this geometry is well defined in the x + > 0, x~ < quadrant. For x + < Xq 



6 -po = e -<Po 



C - X 2 x + x~ -^\og(-X 2 x + x-) , (3.3) 

which is defined for a > a cr , where a = log(— X 2 x + x~) is the static coordinate. At 
cr = a cr there is a singularity, and this line is interpreted as the boundary of the 
strong-coupling region. 

For x + > Xq one obtains an evaporating black hole solution 

M K 

C + — -\ 2 x + {x- + A)--log{-\ 2 x + x-) . (3.4) 
A 4 
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e -Po = e -</>o 



The corresponding Hawking radiation flux at X + is determined in the operator for- 
malism by evaluating 

<* |T__(r)|*o> , (3-5) 

where T (£ - ) is normal ordered with respect to the asymptotically flat coordinates 

^ of the metric (3.4) at X + . Since the £ coordinates are the same as the "out" 
coordinates (2.23) for the classical black hole solution, one obtains by using the point- 
splitting regularization 

r__(r))=^[l-(l + AA e ^)- 2 ] . (3.6) 

The expression (3.6) gives a late-time flux which corresponds to a thermal Hawking 
radiation, with the temperature T = X/2n || |3|. The Hawking radiation shrinks the 
apparent horizon of the solution (3.4), so that the apparent horizon line meets the 
curvature singularity in a finite proper time, at 

*f = ^(-«/4 + eW») • f' ..^ ' < 37 > 

The curvature singularity then becomes naked for x + > xf . However, a static solution 
(3.3) of the form 

e -Po = e -0o = q - \ 2 x + (x- + A) - ~ log(-A 2 x + (x" + A)) (3.8) 

can be continuously matched to (3.4) along x~ = x~ if C = — — log(/c/4)]. 
A small negative energy shock-wave emanates from that point, and for x~ > x^ the 
Hawking radiation stops, while the static geometry (3.8) has a null ADM mass. There 
is again a critical line a = a cr , corresponding to a singularity of the geometry (3.8). 
Note that the scalar curvature of (3.8) is bounded at x~ = x~ , and the singularity 
comes from the pathological behaviour of e~^°, which becomes ill-defined for x~ > x~ . 
This singularity can be interpreted as the boundary of the region where higher-order 
corrections become important. The spatial geometry of the remnant (3.8) is that of 
a semi-infinite throat, extending to the strong coupling region. 

Note that the form of the Hawking flux (3.6) is necessary but not a sufficient 
condition for the Hawking radiation to have a Planckian spectrum. In order to de- 
termine this, one has to evaluate the Bogoliubov coefficients. Since the "in" and 
the "out" coordinates are exactly the same as in the classical black hole background 
case, one can conclude that the one-loop Bogoliubov coefficients are the same as the 
classical Bogoliubov coefficients calculated in ||. Consequently, a Planckian one-loop 
spectrum should be obtained in the late-time approximation. However, by a closer 
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examination of the one-loop geometry one notices that the one-loop horizon is given 
by the line x~ = xj. This means that the classical horizon x~ = —A has undergone 
a small shift due to the one-loop backreaction. The shift —A — x^ is very small, of 
the order of Ae~ M / x . Still, as we are going to show, this shift will have a non-trivial 
consequence for the spectrum of the Hawking radiation. 

The Bogoliubov coefficients can be calculated by using the formalism of || . The 
"in" plain-wave basis is given by 

Uw = * e --T , (3.9) 
V 2lo 

while the "out" basis is given by 

1 



'2w 

where 77, is the position of the one-loop horizon, given by 



e-^flfa-TT) > (3-10) 



Xrji = log 



1 _ ^ e -f (M/A+O-l 

4 



(3.11) 



In the semiclassical approximation, when M >> A, the shift is an exponentially small 
quantity. Still, it produces a change in the classical "out" basis, for which rji = 0. 
Following H, one obtains 

aw = — / dr) Vudrj-ifci , Pu U > = - drj v^-u^ , (3.12) 

7T J— oo 71 J— oo 

so that 

aw = 7^1/- r dr]-e-^-+^- , (3.13) 

while 

1 /a;' rn 



= 2^V/- / ^-e-^"-*""" = -icw . (3.14) 



From (2.23) we obtain 



— r ^- p *flog[AA(e-^ -l)]±iu>' V - 



<£, = f V- / ^"e**^^ ~ 1)J± ^"" , (3.15) 
where a + = a and a - = (3. The integral in (3.15) can be rewritten as 

J 1= r dx(l-x) a x- l - a+b (3.16) 



where a = iu/X, b = ±zu//A, X{ = e Xvi . The integral (3.16) is the incomplete Beta 
function 5(1 + a, 6 — a, sc,) ||10||, so that 



a L' = ^^(^y uj/X B(l+ l co/X,~ l (co T co')/X,x l ) . (3.17) 



The incomplete Beta function can be evaluated for Xi close to one via the expansion 

m 



b-a oo /i i a 

B{l + a,b-a;x) = B{l + a,b-a)-± E / 9 T (1 ~ x T +l+a > ( 3 - 18 ) 

1 + a n=Q (2 + aj n 

where (c) n = When Xj = 1, the formula (3.18) gives the expression for the 

Bogoliubov coefficients obtained in [[J. 

The number of particles of a given energy emmited at 1 + is given by 

POO 

N„ = / Ifl^lW , (3.19) 



which is divergent when (3.17) is used. This divergence is the artefact of the plain- 
wave basis we are using, and can be avoided by using a normalisable basis |4[]. A 
convinient basis is formed by the vawe-packets 

v je {n) = ~r due—^Vu , (3.20) 



which are centered around a~ = ^p, where n,j 6 N |. In the late-time approxim- 
mation, i.e. when n is large, the main contribution in the /^-coefficient integral comes 
from the vicinity of the horizon, for A77 = 0(exp(— 27m/e)), so that one can use the 
approximation e~ Xv — 1 ~ — which gives 

/V - r *T e^-^-^~ . (3.21) 

111 \ UJ J-00 

This is also equivalent to large u' asymptotics of (3.17), since the contribution from 
small rj~ is equivalent to logo/ « 1-nnje. 

The expression (3.21) can be evaluated by using the integral 

(■OO 

Ji= / dxa? a e -ta = 6 _1_a r(l + a, 6e<) , (3.22) 

where T(a,x) is the incomplete gamma-function, a = iu/X, b = ±iu'/X and q = 
—Xf]i. The following expansion 

r(«,x) = r(«)-^^L- , (3.23) 

^0 n!(a; + n) 



can be used to approximate the incomplete gamma function for small x ||10|| , while 
for large values of x one can use [ |10[ 

r(a,x) = x a - 1 e- :E (l + ^^ + 0(x- 2 )) . (3.24) 
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(3.25) 



By using the approximation (3.21) and the formula (3.22) we obtain 

ofc, « ^ X ^(AA)-/ A (±^'/A)- 1 ^ / A r(l + iu/X, Tiu'vi) ■ 

Note that the approximation (3.25) follows from the exact formula (3.17) when oj' is 
large, in accordance with the expected late-time asymptotics. This can be seen from 
(3.18) by approximating each term in (3.18) by its leading large-6 expression. In this 
way one obtains the first two terms of (3.25) when expanded via (3.23). 

From (3.25) and (3.23) we obtain 



a, 



(3, 



r(l + iu/X) -i 



l+icu/X 1 



\l+iu>/\ 



+ 



(3.26) 





2 [ 




~ exp 





r(i + iu/\) + ^(-^,) 1 ^a + . . . 

One can see that the expression (3.26) gives in the limit 77^ — > a Planckian spectrum 
with the temperature T = X/2n. However, the rji corrections do not correspond to a 
Planckian spectrum with a shifted temperature T + 5T, as one would naively expect. 
This can be also seen by applying the large oj' formula (3.24) to (3.25), so that 

where T is independent of oj' . Note that when oj' — > 00 the expression (3.27) tends 
to one, which in the thermal case would mean that T — > 00. This indicates that 
the corresponding Hawking flux will diverge. We can calculate the Hawking flux by 
calculating N u , since 

(T) = r dojojN u . (3.28) 
Jo 

N w can be calculated by going into the basis (3.20), so that 

' ~"""7W , (3.29) 



y/e J j 



U + 2-niv 

due 



where oj 
limit 



je. By applying the expansion (3.23) to /3„ w /, we obtain in the large-n 



|^(n)| 2 = — (e 2 ^ A -l 



_ x 1 sin 2 e2/2 



irXe 



OJ 



+ 



irjie 2A 



7r 2 eoj 



r ( 1 + iuj / X ) e ieti+h(z- 



ioj/X z\z 



c.c. 



+ 0( V 2 ) , (3.30) 

where c.c. denotes the complex conjugate of the previous term, while z = — j log (a// A 2 A) + 
^ and z 1 = \ log(-^A 2 A) + In order to derive (3.30) one needs 

. 2 n.<= ... . Is 

a 
10 



rU 



due mu = - sin — e i{l+ i )ea 



(3.31) 



and 

|r(l + ^)| 2 = -4^ • (3.32) 
sinn 7rx 

The uj' integral over the first term in (3.30) gives the standard thermal contribution. 
However, the uj' integral over the second and higher order in rji terms in (3.30) are 
all divergent. This is caused by the fact that we have taken only a finite number of 
terms in the expansion (3.23), which is a good approximation only for a sufficiently 
small x = u f €i , or equivalently not too big uj' {uj' < e^ 1 ). For big uj' {uj' > ), the 
asymptotic expansion (3.24) is better. Consequently, we have 

poo 

12 , / a, J\ a |2 



N u (n)= I du/|/Uf + J^duj'\l3 u 

POO 

duj'\(3:j 2 + / duj'\PU 2 (3.33) 



o 

where (3 s and (3 l stand for small and large uj' approximations, respectively. The first 
integral in (3.33) has the integrand given by (3.30), and because of the cut-off in uj', 
it is finite. In the second integral we use the approximation (3.24), so that 

1 / 7 —iu>'/\ 

^ = ^V-(AA)^ A ^-^exp(z(^/A)log6 2 ) . (3.34) 

LlX \ UJ —IUJ 

This gives a logarithmic divergent contribution. However, when e« = 0, the large uj' 
asymptotics changes by a phase-factor exp[— i(u/X) log(u//A)], which gives a finite 
N u (n). 

Note that the one-loop flux is also given by the expectation value (3.6), which is 
finite. From the Bogoliubov coefficients analysis it follows that a finite flux is possible 
only if the shift < 0, while the one-loop geometry seems to suggest that e s > 0. 
The way out of this paradox is provided by the fact that the line x~ = x^ gets very 
close to the strong-coupling region. Therefore the one-loop geometry is not a good 
approximation there, and hence our assumption that the effective horizon was given 
by the line x~ = x~ was not a good assumption. 



4. Two-loop spectrum 



A two-loop metric can be obtained by truncating the expansion (2.13) after n = 2. 
Outside of the region occupied by the matter pulse centered around x$, the metric 
can be written as 

ds 2 = —e^ 2 dx + dx~ ; 



,<t>2 



■:,<k> 



1 + e 2 ^ (C_ + C + {x + - x+) 2 6{x + - x+)) 
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(4.1) 



where is the one-loop dilaton solution 

e -*° = C - a(x+ - x£)6(x + - x$) - \ 2 x + x~ - %log\\ 2 x + x-\ , (4.2) 

whith a = X 2 A and C± are integration constants, analogous to the constant C of the 
one-loop solution. The difference now is that C + depends on the matter pulse profile, 
and it is negative for very narrow pulses (those which are shorter than a critical 
length defined by the momentum cut-off ||) and positive otherwise. The constant 
C- is matter independent, and its value is regularization dependent. Consistent 
semiclassical geometries appear for C + > ||. 

As in the one-loop case, the relevant quadrant is x + > 0, x~ < 0. For x + < Xq 
the solution (4.1) is static, and it can be written as 



e 2 ^a 2 



C + e a -ia , (4.3) 



where C_ = —a 2 and a = log(— \ 2 x + x~) is the static coordinate. The solution (4.3) 
describes a two-loop corrected dilaton vacuum. The corresponding scalar curvature 
diverges at the line 

e -*> - a = , (4.4) 

which corresponds to the one-loop singularity line with C replaced by C — a. The 
semiclassical geometry will be defined for a > a c , where C + e ac — f cr c = a. The 
curvature singularity will be absent for C > a — |(1 — log |). Therefore if we want 
to avoid a naked singularity at two loops, the one-loop C has to be increased to 
C = a — |(1 — log |). Note that in the case when C_ = a 2 , there is no curvature 
singularity, and the only singularity comes from the one-loop critical line determined 
by the second formula in (4.3). 

For x + > Xq the solution (4.1) becomes 



2 ^(a 2 -C+{x + -x+) 2 )] , (4.5) 



where 



e 



-00 



C + ^--\ 2 x + (x- +A)-%\og(-\ 2 x + x-) . (4.6) 
A 

It describes a two-loop corrected evaporating black hole geometry. The curvature 
singularity line is given by 

e -20O_ a 2 + C , + (x +_ x + )2 = o _ (47) 



When C + = 0, the equation (4.7) becomes a shock wave singularity equation (4.4), 
and it describes a singularity line of the one-loop metric with a smaller ADM mass 
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C + ¥- — a. The apparent horizon line is given by the equation d + (p2 = 0, which can 
be rewritten as 

x - + A + ^_ = -2C + (*+-s + )e-*VA» (4g) 
4A 2 x+ e- 2 *o - 3a 2 + 3C + (x+ - x+) 2 ' 1 ' } 

In the shock-wave limit (C+ = 0) the intersection point of the apparent horizon and 
the curvature singularity line is given by the one- loop expressions (3.7). The change 
in the intersection point when C + ^ can be evaluated perturbatively |J, and it is 
given by 

6x7= -^ (g<+ ~+° +) W + OOT 
2a axj 

Sxf = - 2X * C+ {xt - X ° )2 x-xf + 0(/3 4 ) . (4.9) 
na a 

The line x~ = Xh = x~ + 5x~ can be considered as the horizon of the two-loop 
semiclassical geometry. Clearly the quantum corrections have shifted the position of 
the horizon, which will affect the Hawking radiation. In the region x~ > Xh a naked 
singularity will appear, unless we impose an appropriate boundary condition. In the 
shock-wave case, one can impose a static solution (4.3) for x~ > x^ 



e <p2 = e <P0 



1 - e*a 2 



, e-^ =C-\ 2 x + (x- + A)-%log{-\ 2 x + {x- + A)) , 

(4.10) 

which can be continuously matched to (4.5) at x~ = Xh if 

£ = ~(l-log£)-a • (4.11) 

However, in contrast to the one- loop case, the metric (4.10) has a curvature singularity 
at a = cr cr , and the naked singularity is not removed. When C + < 0, there is no static 
two-loop dilaton vacuum solution for x~ > Xh which can be continuously matched to 
(4.3) and the naked singularity remains. However, when C + > 0, the naked singularity 
lays in the strong-coupling region, and therefore it can be ignored since it lays outside 
of the region of validity of the metric (4.1). Also note that in the case when C_ > 0, 
the solution with C + > does not have a curvature singularity. As a consequence of 
these properties the Hawking flux will be free of singularities only for C + > 0. 

The asymptotically flat coordinates (a + ,a~) at Zj£ are given by 

\x + = e Xu+ , X(x~ + A) = -e~ xa , (4.12) 

where 

= ■ ( 4 - 13 ) 
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This can be written as 

a- = -ilog[-A(A + x-)] + ^(A + x-)- 2 , (4.14) 

so that a non-zero correction appears at two loops. As a direct consequence, the 
Hawking flux at 1 + will not have a zero-loop form (3.6), which can be seen by eval- 
uating 




48 V rf 

where the primes denote d/da~. One then obtains 

rj- it \\ A 2 j/V (-Ay + 3A 2 ) - C + P a (y) + CjP^y)] 

where y = x~ + A, P 2 = -2y 2 + 3 Ay - 3/2 A 2 and P 4 = y 2 (-^y 2 + 10 Ay - 5A 2 ). 
From (4.16) one can see that T(y) does not diverge for C + > 0. In that case T(y) 
goes to zero for late times (y — > 0), indicating that higher-order loop corrections can 
turn off the Hawking radiation. This also means that the spectrum of the Hawking 
radiation can not be thermal for very late times. 

Relation (4.14) clearly shows that the two- loop "out" coordinates are not the same 
as the analogous zero and one loop coordinates, which indicates that the spectrum 
of the radiation will change. This can be confirmed by evaluating the Bogoliubov 
coefficients in the late-time approximation. The two-loop Bogoliubov coefficients can 
be evaluated from the formulas (3.13-14), where now \a~ = — \og[\A(e~ Xv — 1)] + 

2A 4 A 2 



c^( e - Xr r - l)- 2 and 



Mi = lo § 

so that 



1 _ k-^M/X+O-l 

4 



C + [ % T , (4.17) 

2aax7 



(4.18) 



For late times, we can use the approximation e Xv -1» — A77 , so that 

1 Uj' pH , ( .00 , , , o. _ N lUjCj, 



Ol 1 ss 



1 Uj' fVi f .UJ „ IUJU+ ,_ o , . , \ 

27 V U loo ^ 6XP V A l0g( " A ^ " 2A^ (Arr) ~ ± 77 "J 



(4.19) 



Equation (4.19) can be rewritten as 



' '-(AA) i "/V 2 (a,f.,c) , (4.20) 



2ti\ V 
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where a — iu/X, b — ±iu//A, = — Xr/i and c = iuC + /2X 5 A 2 . Since we are interested 
in the large-6 asymptotics, we can analyze J 2 by using the method of stationary points. 
We can rewrite J 2 as 

POO 

J 2 = / dyy a e-fM > (4.21) 



1 

where /(?/) = by + cy~ 2 . /(?/) has a stationary point at y = (2c/6) 3 , which for a 
sufficiently large b will lay outside of the interval [e;, 00). Hence the main contribution 
to J 2 w iH come from the vicinity of the point y = e*, so that as b — > 00 

j 2 ~-k e *) aex p(-^- ce r 2 ) • ( 4 - 22 ) 



Note that the argument we have used to derive (4.22) clearly holds for a, b, c real 
and positive. When a, b, c are complex, an analogous argument applies (see JTTJ and 



the appendix), and (4.22) is still valid. We also show in the appendix that the same 
method can be applied to the exact expression (4.18), and the result (4.22) is again 
obtained when b is large. 

The expression (4.22) has the same asymptotics as the incomplete gamma func- 
tion (3.22), and therefore we will encounter the same problem as in the one-loop case. 
Since N u (n) will then diverge, the corresponding Hawking flux will also diverge for 
late times. On the other hand, the Hawking flux from the expectation value (4.16) is 
finite for late times, and we have the same discrepancy as in the one-loop case. Since 
the line x~ = xC + 8xi~, which is suggested by the semiclassical geometry as the 
horizon, gets close to the border of the strong-coupling region where the spacetime 
geometry is not well-defined, we can invoke the same argument as in the one-loop 
case, and resolve the discrepancy by taking < 0. In that case the stationary-point 
method gives thermal Bogoliubov coefficients (see the appendix), and we obtain a 
Planckian late-time spectrum in the leading-order approximation with the temper- 
ature T = 3A/27T. This increase in the temperature is an expected back- reaction 
effect. However, an unexpected feature is that the increase is not perturbative in C + . 
A related problem is that the Hawking flux (4.16) deviates only slightly for small 
enough C + from the zero-loop flux which corresponds to the temperature T = X/2tt. 
Note that for big enough C + there is a late time y§ for which T(y Q ) coincides with 
the thermal flux corresponding to the temperature T = 3X/2tt. In this case one can 
argue that there is an agreement with the Bogoliubov coefficients calculation when 
y ~ yo, while for the times y > y , the discrepancy can be justified by the argument 
that the background geometry is not well-defined (strong-coupling region) and hence 
the Bogoliubov coefficients are not well defined there. Still, there is a discrepancy in 
the weak-coupling region for small C+, and a possible explanation is that the approx- 
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imation (A. 15) is not good enough, and that the a% and higher-order terms should 
have been taken into account. 

5. Conclusions 

The results of our analysis indicate that small back-reaction corrections in the 
background geometry can induce large changes in the Hawking radiation spectrum. 
We have seen that a small positive shift in the position of the classical horizon gives 
a non-Planckian late-time Hawking radiation spectrum whose flux diverges. On the 
other hand, a small negative or a null shift gives a Planckian late-time spectrum. 
Consistency considerations then select a null or a negative shift as an effective horizon 
shift, contrary to what is suggested by the semiclassical geometry. This is explained 
by the breakdown of the semiclassical approximation for late times near the classical 
horizon. 

In the one-loop case the Bogoliubov coefficients give a late-time Planckian spec- 
trum which is consistent with the operator quantization Hawking flux only for times 
x~ < . After that the BPP solution gives a zero flux, which indicates that one 
enters a non-thermal regime. It is plausibile that the higher-order stationary point 
contributions to the Bogoliubov coefficients may account for this. In the two-loop 
case the relation between the late-time spectrum and the operator quantization flux 
is similar but slightly more complex. The consistency again requires that the time is 
not too late, and the new requirement is that the constant C + is not too small. It 
is likely that this restriction is caused by the low-order approximation we are using 
to calculate the two-loop Bogoliubov coefficients. Including the higher-order terms 
in the stationary-point approximation may also account for a non-Planckian nature 
of the very late time spectrum, which is suggested by the behaviour of the operator 
quantization Hawking flux T . However, a further investigation is necessary in order 
to clarify these issues. 

Note that when C + < 0, there is a total disagreement between the flux from the 
Bogoliubov coefficients and the operator quantization flux T. In that case T diverges 
to minus infinity for late times. This cannot be reconciled neither with the positive 
horizon shift case (the corresponding flux diverges to plus infinity), nor with the 
negative horizon shift case (the corresponding flux is finite and positive). However, 
C + < solution can be considered as an unphysical solution, since it requires a matter 
pulse whose width is shorter than the 2d analog of the Planck length II . 
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Appendix 



In this appendix we give some relevant theorems from |Tl[ , and use them to obtain 
the asymptotic expressions for the relevant integrals. 

In the one-loop case, the relevant integral is 

I^b) = dTT b - a -\l-T) a (Al) 



where Xi = exp(— e,), a = iu/X, b = This integral can be rewritten as 

I 1 (b)=x\- a dte^ie 1 -x.Y (A2) 



o 



where r = x^e *. Now we use a theorem from (TTJ, page 108 

Theorem (1): Let I(z) = / °° dte~ zt q(t), such that q(t) is real or complex, < 
A s e a ^ (a E R) and q(t) holomorphic in the sector S = {ax < Argt < a 2 } (at < 0, 
«2 > 0). Then as z — > oo 

in the sector S n = {— a\ — ir/2 < Arg z < —a 2 + tt/2}. 

In the case of the integral (A. 2), q(t) = (e* — Xi) a , and the conditions of the 
theorem (1) are satisfied, so that as b — > oo 

h{b)^x\- a {l ~* l)a +0(b- 2 ) . (A4) 

Note that the expression (A. 4) has the same large-6 asymptotics as (3.22) for Xj close 
to one. 



When Xi = 1, then a generalization of the theorem (1) applies |TT[|, page 114 
Theorem (2): Let q(t) be holomorphic in S and q(t) = 0(e CT '*') as t — > oo. Let 

oo 

q(t) « a s t {s+1/ ~' x)/ ' x 
for t —>■ + and t G S, where \x > and Re v > 0. Then 



Er — -Ar (^4-5) 
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as z — > oo in the sector SV- 



In our case, q(t) = (e* — l) a « E^Lo a st s+a , so that // = 1 and z/ = l + a = l + zcj/A, 
and hence we can use the theorem (2). Therefore, as b — > oo 



(1 + 0(6-) 

whose leading-order term is (3.22) for q = 0. 
In the two-loop case, we have 



(A6) 



h(b) 



d 



T T 



_6-o-l, 



1 -r) a exp [-c(r -1 - 1)" 



6— a 



dte (e — Xi) a exp 



-cfx7 1 e* - 1 



(A7) 



When Xi < 1, then 



= (e* -Xi) a exp 



-c(x7V - 1)~ 2 



so that the theorem (1) applies, and hence 



— c 



1 - a;, 



s=0 



+ 0{b~ 



(Ai 



It is easy to see that the large-6 asymptotics of the expressions (A. 8) and (4.22) 
coincide for Xi close to one. 

When Xi = 1, then q(t) ~ t a exp(— ct~ 2 ) as t — > + , so that neither theorem (1) 
nor theorem (2) can be used. One can then use a theorem from |IT|, pg 127 

Theorem (3): Let 

I(z)= f dte- zp{t) q{t) = [ dte- zp{t) q(t) 
Ja JC ab 

such that 

(i) p(t) and q(t) are holomorphic and single-valued in a domain T 

(ii) curve C a b G T 

(hi) p'(t) has a simple zero t £ Ca6 

(iv) 6* = Ar^z G [#i,#2L M > ^, $2 — #i < 7T, Z > 0, converges absolutely 
and uniformly at a and b. 

(v) Re{e ie (p(t) — p(to))} is positive on C a fe, except at to, and is bounded away 
from zero uniformly with respect to 9 as t — > a or 6 along C a b- 
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Then as z — > oo 



«2s 



7(z)«2e-^)^r( S + l/2)- , 
for 6*i < Agz < 2) where 



(A9) 



a 



a 2 



„ 2p"V /5p'" 2 p iw 



,7 9 



(2p 



/A -3/2 



J9" \ 6p" 2 2j9' 

and so on. 

In the case of the integral (A. 7) when Xi = 1 one can take 



p(*)=f + -(e'-l)- 2 , 5(t) = (e«-l) a . 



(A10) 



(All) 



However, theorem (3) cannot be applied directly, since the role of z is played by b, 
and in the choice (A. 11) p(t) depends on b. In order to obtain a proper p(t), note 
that if we rewrite J 2 as 

POO 

UK) = / dte~ m (A 12) 

Jo 

where f(t) = — a\og(e t — 1) + bt + c(e* — 1)~ 2 , then a main contribution to I 2 comes 
from the stationary point /'(to) = 0, which for large b is given as to ~ (2c/6) 1//3 . Since 
then t is close to zero, one can approximate f(t) with its small t asymptotics, which 
is the late-time approximation we used in sections 3 and 4. Therefore as b — > oo 



bp(t) ^bt + ct~ l , q(t) « t a , 
l 

and we make a rescaling t = r(c/6) 3 , so that 

6p(t) w zp{r) = z(t + r~ 2 ) , 

where z = (& 2 c) 1//3 . Consequently 

h « (c/6) (a+1)/3 /°° rfr e -^ T+T " 2 V a , 
•/ o 



(A13) 



(A14) 



and the theorem (3) can be now applied. We then obtain in the leading-order ap- 
proximation (ao term) 



/I 



-l) a / 3 | 2 = exp(2cu7r/3A) , 



(A15) 



so that the leading-order two- loop Hawking temperature is given by T = 3A/27T. 
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